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We study superconductivity in spin-orbit-coupled systems in the vicinity of inversion symmetry 
breaking. We find that due to the presence of spin-orbit coupling, fluctuations of the incipient parity¬ 
breaking order generate an attractive pairing interaction in an odd-parity pairing channel, which 
competes with the s-wave pairing. We show that applying a Zeeman field suppresses the s-wave 
pairing and promotes the odd-parity superconducting state. Our work provides a new mechanism 
for odd-parity pairing and opens a route to novel topological superconductivity. 

PACS numbers: 74.20.Rp, 74.20.Mn, 74.45.+C 


Over the last few years, the search for unconventional 
superconductors has received a new impetus from the 
study of topological phases of matter. After early works 
on superfluid Helium-3 [I:] and recent developments on 
topological insulators, it has been theoretically estab¬ 
lished [2| that superconducting states can be classified by 
their topological properties. Unlike conventional s-wave 
superconductors, topological superconductors are pre¬ 
dicted to harbor exotic quasiparticle excitations on the 
boundary. There is currently intensive effort searching 
for topological superconductivity in naturally-occurring 
materials (see for example, Ref. [3 h 13[), though definitive 
experimental evidence is still lacking. For the majority 
of superconductors that are time-reversal and inversion 
symmetric, the single most important requirement for 
being topologically nontrivial is that the pairin g o rder 
parameter must be odd under spatial inversion [3, |l4[ , 
e.g., p- or /-wave. This rekindles interest in finding odd- 
parity superconductors. 

The parity of the pairing order parameter is tied with 
its spin. In the absence of spin-orbit coupling, odd-parity 
pairing is spin-triplet and vice versa. It has long been 
known from studies on superfluid He [151 ] that triplet 
pairing can be induced by enhanced ferromagnetic spin 
fluctuations in the vicinity of ferromagnetic instability. 
This mechanism for triplet pairing, if realized in the solid- 
state, can lead to a topological superconductor 0 , anal¬ 
ogous to the topological superfluid He. 


In this Letter, we propose an alternative mecha¬ 
nism for odd-parity pairing in the vicinity of nonmag¬ 
netic, inversion-symmetry-breaking phases in spin-orbit- 
coupled systems. In the presence of spin-orbit interac¬ 
tion, such parity-breaking orders directly couple to elec¬ 
tron’s spin texture on the Fermi surface 0- As a re¬ 
sult, the fluctuations of an incipient parity-breaking or¬ 
der, which we call “parity fluctuations”, generate an ef¬ 
fective interaction that is strongly momentum and spin 
dependent. Without assuming any special features of the 
Fermi surface, we show on general ground that this ef¬ 
fective interaction is attractive in both the s-wave and 
an odd-parity pairing channel. Moreover, the pairing in¬ 


teractions in the two channels are found to be of the 
same order of magnitude, and in several cases, roughly 
equal. We show that either Coulomb interaction or Zee- 
man field suppresses the s-wave pairing and promotes the 
odd-parity superconducting state on the border of parity¬ 
breaking order. Finally, we propose the pyrochlore oxide 
Cd 2 Re 2 C >7 and doped SrTiC >3 heterostructures as candi¬ 
date systems where odd-parity superconductivity medi¬ 
tated by parity fluctuations may be realized. 

In this work, we consider parity-breaking orders that 
are time-reversal invariant and carry zero momentum. 
Such order may originate from an unstable odd-parity 
phonon or the electron-electron interaction. The order 
parameter can be represented by a Hermitian fermion 
bilinear operator Q with the same symmetry, which takes 
the form 

Q= F^(k)4 a c k/3 , with r f (k) = r(k). ( 1 ) 

k,a/3 


Here a , /3 are pseudospin indices for the two degenerate 
states at every k. In spin-orbit-coupled systems, these 
states are not spin eigenstates, but remain degenerate 
in the presence of time-reversal (0) and inversion (P) 
symmetry fl8| . For simplicity of notation, we have chosen 
an Ising-type parity-breaking order in Eq. ©■ Vector and 
high-rank tensor orders are described by a multiplet of 
Hermitian operators denoted by ; these orders will be 
encountered later. 

Different types of parity-breaking orders are classified 
by their transformation properties under crystal symme¬ 
try operations, which act on electron’s spin and momen¬ 
tum jointly. Before proceeding to the symmetry analysis, 
we emphasize that the form factors for physical observ¬ 
ables, such as T(k) for Q, depend on the basis for the 
doubly degenerate energy band. For the purpose of this 
work, it is most convenient to choose the “manifestly co¬ 
variant Bloch basis” (MCBB) 0- In this basis, the two- 
component electron operator (c^-,, c)^, 2 ) transforms simply 
as a spinor field in k space under time reversal and crystal 
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symmetry operation g £ 0(3): 

® C ka® 1 = £ a/3 C -k/3 (^) 

gcLs -1 = U a p(z )cj c , /3 , (3) 

where k* = gk is the star of k, e a p is the Levi-Civita 
symbol, and U( g) is the U(2) matrix that represents the 
action of g on the pseudospin, in the same way as it acts 
on the spin of a free electron. 

It then follows from the symmetry transformation laws 
mm that the form factor T(k) of time-reversal-invariant 
and parity-breaking orders satisfies the condition 

r(k) = er*(k)e = -r(-k), (4) 

and hence takes the general form E3 

T(k) = dk ■ er, with dk = — d k (5) 

where cr = (a x ,a v ,a z ) denotes Pauli matrices in pseu¬ 
dospin space. The d -vector field dk defines the pseu¬ 
dospin splitting in the ordered state, whose magnitude 
and direction vary over the Fermi surface. 

It may seem counter-intuitive that nonmagnetic parity¬ 
breaking orders, such as structural distortion and or¬ 
bital order, couple to electron’s spin. As we show by 
example in Supplementary Material jlii |. this remark¬ 
able fact is a general consequence of spin-orbit interac¬ 
tion in centrosymmetric systems. It will play a crucial 
role in meditating superconductivity in the vicinity of 
parity-breaking order. In contrast, for spin-rotationally- 
invariant systems, the above symmetry analysis implies 


that parity-breaking orders at zero momentum cannot 
couple directly to electrons on the Fermi surface, unlike 
the nematic order that is even-parity USUI. There¬ 
fore, spin-orbit coupling is crucial for superconductivity 
meditated by odd-parity phonons or parity fluctuations. 

In a system close to a parity-breaking instability, the 
effective interaction arising from the order parameter 
fluctuations is given by 

Hen = ^2 ( 6 ) 

q 

where Q(q) = Q) is the Fourier transform of the 

order parameter field in real space: 

<2(q) = \ ( r «/3( k + q) + r a/ 3(k))4 +qQ c k ,8 (7) 

k ,a/3 

Within the random-phase approximation (RPA), V q can 
be expressed in terms of the q dependent susceptibility: 
Vq = 1/(1 + x(q)/). Vq is enhanced and has a maximum 
at q = 0 close to a q = 0 instability. Restricting the ef¬ 
fective interaction © to the Cooper pairing channel with 
zero total momentum, we obtain the pairing interaction 

Hp — ^2 Vafins (k, k )Ck a C_j c ^gC_k / 'yCk '<5* (8) 

k,k' 

Using (0. © and 0, we find the momentum- and 
pseudospin-dependent interaction vertex V a/ 3 -ys(k, k') is 
given by 


I4/3 7 ,5(k, k') = — - ^Vk_k'(di 


k T dk') * cr a s (dk T dk') * (Tg 7 Vk+k'(dk dk') * rr a7 (dk dk') * crps'j • 

I- 


( 9 ) 


To proceed, we expand Vk±k' hi the pairing interaction 
© in terms of spherical harmonics on the Fermi surface: 
Vk±k' = Vo T Uik • k' + .... Below we consider the lead¬ 
ing term Vo- Despite Vo is a constant, the interaction 
vertex © inherits the form factor of the parity-breaking 
order parameter Q , which is strongly pseudospin- and 
momentum-dependent. It consists of two types of terms: 

V = V e + V°, where V e contains the product of compo¬ 
nents with the same momentum: 

V^sOqk') = -y J2( d l d k + d k'd J k ')KuJ 7 - cr a 7 cr ^)’ 

and V° contains the cross terms: 

v;v(k.k') = -y £(44' + 44)K^ 7 + «). 

hi 


Note that V e (V°) is an even (odd) function of k, k', and 
antisymmetric (symmetric) under exchanging the pseu¬ 
dospin indices either a/3 or qd. Therefore, V e and V° cor¬ 
respond to the even-parity pseudospin-singlet and odd- 
parity pseudospin-triplet pairing channels, respectively. 

The above pairing interaction can be decomposed into 
different superconducting channels that belong to differ¬ 
ent representations of the crystal symmetry group. Be¬ 
fore proceeding, we describe the general classification 
of time-reversal-invariant superconducting order param¬ 
eters, taking the form 

^ = 2 ^2 e / 3 T-f“^( k ) C ka C -k 7 ' (19) 

k,Q' y 07 

where the form factor F(k) satisfies the symmetry con- 
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dition 


F*(k) =F{k) = -eF*{- k)c. 


( 11 ) 


Moreover, it follows from ([3]) that e^ 7 cj ca cL k7 has the 
same transformation law under crystal symmetry op¬ 
erations as c kQ Ck^- This implies every time-reversal- 
invariant superconducting order parameter has a coun¬ 
terpart in the particle-hole channel, with the same sym¬ 
metry. In particular, odd-parity superconducting order 
parameters, which have F(k) = — F(— k), admit the 
same classification as particle-hole orders Q described 
earlier. 

To proceed with the classification, it is instructive to 
first consider the most symmetric group 0(3), the group 
of all joint 3D rotations and reflections of spin and mo¬ 
mentum, from which all point groups descend. In this 
case, all possible odd-parity orders defined by the form 
factor © are classified by the total angular momentum 
J and the orbital angular momentum L (which must be 
odd) E3- At the lowest order L = 1, there are three 
types of particle-hole orders: gyrotropic, ferroelectric and 
multipolar. The corresponding form factors are listed in 
Table I. Classification for 2D systems with 0(2) symme¬ 
try is also presented. 

As expected from the one-to-one correspondence be¬ 
tween particle-hole and particle-particle orders, these 
form factors also classify odd-parity pairing symmetries 
of spin-orbit-coupled superconductors. For example, the 
pairing order parameter with the isotropic form factor 
Ti = k cr coincides with a particular choice of order 
parameters for the Balian-Werthamer phase of He-3 [22| . 
On the other hand, the pairing order parameters with 
anisotropic form factors T 2 and T 3 are time-reversal in¬ 
variant and spontaneously break the rotational symme¬ 
try, resulting in an odd-parity superconductor with ne¬ 
matic order [23} • To our knowledge, such anisotropic 
phases have not been found in He-3; their existence re¬ 
quires spin-orbit coupling. 

We now use the effective interaction H e g given by 
Eq.©) to study superconductivity in the vicinity of 
each type of parity-breaking order in Table I; for multi- 
component operator summation over fj, is taken. In 
all cases, we restrict H e g into Cooper pairing channel 
with zero momentum, and decompose the pairing inter¬ 
action H p into various superconducting channels: 


Hp = V 0 (a 0 S ] S + E“"E^)» 


( 12 ) 


where & = (1/2) J2k, a p e a/ 3 c k,a c - k/3 is the s-wave su¬ 
perconducting order parameter, and F£' denotes various 
odd-parity order parameters defined in Eq. (flOl) and clas¬ 
sified in Table I. Coefficients a n take different values for 
different types of interactions, and all are gathered in the 
Table El Details of our calculation can be found in Sup¬ 
plementary Material 0 - Since Vq < 0, a n > 0 means 


3D system with 0(3) symmetry 

transformation property 

Tr(k) = (k-<r) 

pseudoscalar 

It(k) = [k x or]* 

vector 

F^(k) = fcV + fcV - i(k ■ <r)6 ij 

rank 2 tensor 


2D system with 0(2) symmetry 

transformation property 

Ti(k) = k x a x + k y a y 

pseudoscalar 

r 2 (k) = k x a y - k y a x 

pseudoscalar 

ri(k) = 

vector 

r:/: k) = fcV + fcV - (k • cr)5 ij 

rank 2 tensor 


TABLE I: Classification of odd-parity order parameters for 
spin-orbit coupled systems, and their transformation proper¬ 
ties under joint spin and momentum rotations in three and 
two dimensions. Rank 2 tensor 1^3 (IA) is symmetric and 
traceless, and hence has 5 (2) independent components. 


type of interaction 

a o 

ai 

02 

03 

Qi(q)Qi(q) 

1 

1 

-1 

0 

Oa(q)QS(-q) 

2 

-4/3 

1/2 

-1/4 

Q!/(q)Q(3 l (-q) 

20/3 

0 

-5/3 

-1/2 


type of interaction 

Oo 

Ol 

02 

03 

a4 

Qi(q)Qi(q) 

1 

1 

-1 

-1 

0 

Q2(q)Q2(-q) 

1 

-1 

1 

-1 

0 

Qs(q)Q l 3(-q) 

1 

-1/2 

-1/2 

1 

- 1/4 

Q 4 (q)Qr(-q) 

4 

0 

0 

-4 

0 


TABLE II: Decomposition of the different types of interaction 
into different pairing channels, see Eq. m- n = 0 denotes 
the s-wave channel; n = 1,..., 4 denotes the odd-parity chan¬ 
nels classified in Table I. a„ > 0 corresponds to attractive 
pairing interaction. 


attractive interaction in the corresponding pairing chan¬ 
nel. 

From Table El we obtain the superconducting insta¬ 
bility driven by each type of parity fluctuations. In all 
cases, there is an instability in the s-wave channel, sim¬ 
ilar to phonon-meditated pairing in conventional super¬ 
conductors. More importantly, in all cases except the 
multipolar orders T 3 and T 4 , there is also an instabil¬ 
ity in the odd-parity channel with the same symmetry as 
the incipient particle-hole order that drives superconduc¬ 
tivity. Remarkably, for Ising type orders described by a 
single-component Q, the pairing attraction in the odd- 
parity channel is equal to the one in the s-wave channel, 
leading to identical superconducting transition tempera¬ 
tures. For ferroelectric type orders described by a vector 
Q l , the pairing attraction in the odd-parity channel is 
weaker than, but still of the same order of magnitude as, 
the one in the s-wave channel. 

Although fluctuations of multipolar orders in rotation- 
ally invariant systems do not lead to pairing in any odd- 
parity channel, the situation becomes different in real ma- 
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terials where the crystal symmetry is taken into account. 
In any crystals, the five components of rank 2 tensor 
invariably split into more than one representations of the 
point group. For example, for Oh point group, the diag¬ 
onal and off-diagonal components of Ts split to form e g 
and t 2 g representations, which have 2 and 3 independent 
components respectively. For many point groups such as 
D 4 / 1 , T 4 also splits into one-dimensional representations, 
with form factors k x a x — k y a y and k x a y + k y a x respec¬ 
tively. We find fluctuations of such multipolar orders of 
reduced symmetry generate attractive pairing interaction 
in the odd-parity channel of the same symmetry. The in¬ 
teraction strength is weaker than the s-wave channel in 
the case of e g and t 2g orders, and is equal to the latter 
in the case of Ising type f 4 orders [h|. 

The above finding of odd-parity pairing meditated by 
parity fluctuations in spin-orbit-coupled systems is the 
main result of this work. It is interesting to make a com¬ 
parison with the mechanism of triplet pairing meditated 
by spin fluctuations. In that case, the effective interac¬ 
tion is given by V (q)s(q) • s(— q), where s is the spin 

operator and V(q) = 1/(1 + I\ s ( q)) is determined by 
the spin susceptibility y s (q). Importantly, to obtain the 
pairing interaction in the triplet channel requires (q) to 
have a nontrivial q-dependence. Approximating \ s (q) by 
its zeroth spherical harmonic, which is a constant, does 
not generate triplet pairing, simply because two electrons 
at the same spatial location cannot form a triplet. In 
contrast, we obtained odd-parity pairing in this leading- 
order approximation, without relying on any special fea¬ 
tures of the susceptibility of parity-breaking order. 


Given that the pairing interaction we found has compa¬ 
rable or even identical strengths in the s-wave and odd- 
parity channels, small residual interactions or external 
perturbations become important in lifting the degener¬ 
acy and eventually determine which one of the two com¬ 
peting pairing symmetries is realized. An in-depth study 
of the effects of residual interactions necessarily involve 
material-specific details, which is beyond the scope of 
this work. Nonetheless, it should be noted that Coulomb 
interaction is most repulsive and pair-breaking in the s- 
wave channel, which can make the odd-parity pairing en¬ 
ergetically favorable. This role of Coulomb interaction in 
the competition between s-wave and odd-parity pairings 
has been recognized [ 3 } and emphasized 241 in recent 
model studies. 


In addition to Coulomb interaction, the s-wave pair¬ 
ing is suppressed by a magnetic field B that splits the 
spin degeneracy, which is pair-breaking and sets the Pauli 
limit for the upper critical field. However, Zeeman spin 
splitting has variable effects on odd-parity superconduct¬ 
ing states in spin-orbit-coupled systems, as we show now. 
First, let us consider how the doubly degenerate bands at 
every k, or the pseudospin, split under a Zeeman field. 
The coupling of pseudospin to Zeeman field takes the 


general form 


H z = Y c] ! _gij(k.)B i a j (k)c k . 


(13) 


The g- factor gi g ( k) is a function of k, and can be ex¬ 
panded into different spherical harmonics over the Fermi 
surface. Importantly, since the pseudospin operator ay is 
defined in the manifestly covariant Bloch basis and has 
the same symmetry as electron’s spin, gij( k) generally 
has a dominant zeroth spherical harmonic component 
g~j . Assuming gij (k) ~ g/-. we obtain a uniform spin 
splitting over the Fermi surface, with a spin quantization 
axis in the direction of hi = g/jBi. The Pauli limit will 
be absent for the odd-parity pairing if its d-vector d(k) 
is perpendicular to h, for all k on the Fermi surface. 
For example, in 2D systems with rotational symmetry, 
an in-plane field B induces a spin splitting in the direc¬ 
tion parallel to the field. The odd-parity pairing with 
T 3 (k) = (k x a z ,k y a z ), whose d-vector is out of plane, is 
not Pauli limited. Therefore, Zeeman field is an effective 
way of tuning the competition between different pairing 
symmetries and promoting certain types of odd-parity 
superconductivity for which the Pauli limit is absent or 
largely enhanced. 

Finally, we propose candidate materials for odd-parity 
superconductivity in the vicinity of parity-breaking or¬ 
der. First, the pyrochlore oxide Cd 2 Re 2 07 undergoes 
a continuous parity-breaking phase transition at T p = 
200K (2f| [26| with a large mass enhancement of con¬ 
duction electrons |27|, and becomes superconducting at 
T c = 1.1K [28|. The application of high pressure has 
significant effects on these phases, and generates a vari¬ 
ety of new phases identified from resistivity anomalies. 
Remarkably, around a critical pressure of P c = 4.2GPa 
where the parity-breaking order is suppressed, an anoma¬ 
lously large upper critical field of 7.8T is observed, which 
is 27 times larger than at ambient pressure and signif¬ 
icantly higher than the Pauli limit 4.2T evaluated as 
H p = 1.84T C [2§| . These phenomena seem to fit into 
the theoretical picture presented in this work. Therefore, 
we propose that the superconducting state of Cd 2 Re 207 
around P c is driven by parity fluctuations, and may have 
an odd-parity pairing symmetry. 

Another candidate system is inversion-symmetric het¬ 
erostructure of doped SrTiOs with intrinsic spin-orbit 
coupling [36j. Bulk SrTiOs is close to the ferroelectric 
instability and becomes superconducting upon electron 
doping [j37, 38|. In doped SrTiOs heterostructures with 
the superconducting dopant layer of a few nanometers 
thickness, the in-plane upper critical field exceeds the 
conventional Pauli limit [36|. It is worthwhile to ex¬ 
amine the possibility of an odd-parity superconducting 
state under a large in-plane field, as we discussed ear¬ 
lier. A model study for possible superconducting phases 
in SrTiOs heterostructure will be presented elsewhere. 
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Throughout this work, we have stayed away from the 
immediate neighborhood of the quantum phase tran¬ 
sition point, where long-wavelength and low-frequency 
fluctuations of the parity-breaking order pile up and the 
RPA type effective interaction used in this work is inap¬ 
plicable. The physics in the quantum critical regime is 
an interesting topic which is left to future study. 
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I 

Supplemental materials 

This Supplementary Material consists of two sections. In section I, within an inversion-symmetric two-orbital model, 
we explicitly derive the Fermi surface form factor of a parity-breaking order, which elucidates the important role of 
spin-orbit coupling. In section II, we provide a detailed derivation of the decomposition of effective interactions into 
different pairing channels. 


I. BILAYER RASHBA MODEL 


We consider a two-layer system with intrinsic spin-orbit coupling, which is the 2D analog of the model considered by 
Fu and Berg for the topological insulator Bi 2 Se 3 [Sll j. The same model applies to inversion-symmetric heterostructures 


of SrTiC >3 [S2|. Importantly, the two layers in the system considered here, denoted by t z below, are interchanged 


under the inversion. It follows from symmetry that electrons on a given plane experience an out-of-plane electric field, 
which is opposite for the upper and lower plane. This local electric field generates local Rashba spin-orbit coupling 
for electron’s motion within the plane, with opposite signs for the upper and lower planes. In addition, there is an 
inter-plane tunneling, which takes the off-diagonal form mT x . Based on symmetry considerations, the Hamiltonian 
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2 -band model 

1 -band model 

t v s z 

Fi(k) = k x a x + k y a y 

T Z 

F 2 { k) = k x a y - k v a x 

T y S X 

Ef(k) = k x a z 

T v s v 

F y {k) = k y a z 


TABLE SI: The correspondence between pairing orders in 2-band model considered in fsTIl and 1-band model presented in this 
Letter. 


for such a bilayer Rashba system, up to first order in k, is given by: 

r (l\r 

H = J (27t ) 2 [ m7 ~ X + v(k x s v - k v s x )t z - jn]V> k- (SI) 

where p is chemical potential, r* and s l are Pauli matrices in the layer and spin space, correspondingly. Without loss 
of generality, we assume the chemical potential lies in the conduction band. 

Since the superconducting gap is much less than chemical potential, it is convenient to consider the conduction 
band only, and completely ignore valence band. To do that, we need to expand creation and annihilation operators 
V’k ; V’k hr terms of the eigenstates of the conduction and valence bands, and then simply omit the contribution of the 
valence band. This expansion takes the form 


V’k = 5Z a ka C ka: 

ori 


(S2) 


where i = c,v labels conduction/ valence band, and a = 1,2 denote band eigenstates in the MCBB. In this ba¬ 
sis, eigenvectors a£ Q corresponding to the states of the conduction band, which is the two-orbital analog of Bloch 
wavefunctions in the continuum, have the form 


( A ^ 


/ -ik-P- \ 

ik+j3- 


P+ 

P+ 

) a k2 ~ 

ik-P- 

V -ik+P- y 


V P+ / 


(S3) 


where k± = k x ± ik y , and /3± = (l/2)y/l ± (m//z), /r = \Jm 2 + v 2 k 2 F . Mathematically, the mapping onto conducting 
band simply implies that in (IS2I) we keep a£ a only, and omit terms with a^ a . 

Now, we consider an inversion-symmetry-breaking charge order: an electron density imbalance on the two layers, 
(ni — ri 2 )(r) = i/> + (r)r z if>(r) , which does not involve electron’s spin. After the mapping onto the conduction band, 
this expression can be rewritten in the form 

V’ + (r)r^(r) -A ^ e ^' C ka(^+P) V a/3 C P/3 e_ir ' (k “ P) > ( S4 ) 

k.p 


which exactly coincides with the Fourier transform of T 2 (see Table H]). Now, we see explicitly this form factor 
is pseudo-spin dependent. This derivation illustrates how inversion-symmetry-breaking charge order in spin-orbit- 
coupled systems couples to the pseudospin degree of freedom on the Fermi surface. 

It is also instructive to present the correspondence between p-wave pairing order parameters in single-band model 


and in 2-band model considered in jSll . lS3l . lS7l ]. We demonstrate the derivation of the correspondence on the example of 


the two-component order parameter, A^ = V’ka(*'r y s I s y ) Q ,/3'i/'_k/3. To establish the correspondence we, again, perform 
mapping onto conduction band using MCBB, Eq. (IS3l) . The result of this mapping reads as follows: 


Afc = i’k a {iT V S l S V ) a p'lp^ ± -k l [ck\C-k2 + C k 2C-kl] = ±—^-Cka(*fcVV y ) a/ 3C_ k /3, ( S5 ) 

M A 4 

where + (—) sign corresponds to i = x (i = y). The correspondence between other pairing orders can be obtained 
analogously. The result is shown in Table [ST] 
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II. DECOMPOSITION OF THE EFFECTIVE INTERACTION INTO SUPERCONDUCTING CHANNELS 

Here we show the detailed derivation of Eq. m and Table EH We start with the effective interaction f9]) and use 
the approximation V(k ± k') = Vo- We consider all possible types of interactions, i.e. gyrotropic, ferroelectric and 
multipolar, in both three and two dimensions. 


A. Rotationally invariant interaction 

We start with the case of rotationally invariant interaction. We demonstrate the details of the calculation using 
gyrotropic interaction in 3D as an example, and then present the results for all other types of interaction in 3D and 
2D. 

The task is to decompose the effective interaction into different p-wave superconducting channels, classified in Table 
El and, possibly, s-wave channel. Mathematically, it means that interaction matrix elements given by Eq. |9]) need to 
be rewritten in the form: 


V afh a{ k,k') = E (a 0 (ia"Mi<7% + a 1 (iF 1 (kK)^(*F 1 (k>% + a 2 {iF' 2 (K)a v ) a p{iF^)cj y )+ s + 

+a 3 (iF^ J (k)a v ) a0 {iFi\k')a y )+ s ^ , (S6) 

where summation over repeated indices i, j = x,y,z is implied. All F,; have the same momentum and spin dependence 
as Ti in Table El but denote pairing in particle-particle channels. 

For the case of gyrotropic interaction, the Hamiltonian has the form 


HeS — 4 a ka a k'-q/3 ak '7 a k-q<5 ^ 2 k • CT a $ ^ 2 k ■ (T. 


(S7) 


k.k'.q 


Restricting this Hamiltonian to the Cooper channel only, we obtain the pairing Hamiltonian (j8j) with the vertex 


V a p l5 (k, k') = ((k + k') • cr a5 (k + k') • cr^ - (k - k') • cr a7 (k - k') • crps^j . (S8) 

This interaction corresponds to the choice = k in Eq. ©■ 

As we pointed out earlier, this vertex consists of two different terms. The first one contains the product of 
components of the same momenta: 


^(k, k') = (S9) 

ij 

This term describes s-wave pairing and corresponds to the first term in Eq. (|S6[) . Indeed, using the identity 

<6~ ~ vfaKy = -2(if7 s/ ) a ^(icr 2/ )+ (5 ^ (S10) 

we find 


y a e j97i (k,k / ) = ^(*^) aj8 (*^)+, 

therefore, ao = 1, and we have attraction in the s-wave pairing channel. 

Next, the second term in the vertex is the bilinear function of different momenta and has the form 


(Sll) 


= -^ £(***" + + aly 05 ). 


Vfasi k,k') 


(S12) 
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This term is responsible for the p-wave pairing, and need to be decomposed into p-wave superconducting channels, 
listed in Table [I] After some algebra, we find 

7 «(k, k') = y |(i(k • cr)a v ) a p(i(k' ■ cr)a y )+ g - (*[k x <r\a v ) a p(i[. k' x cr]<j y )+ s }. (S13) 

After summation, V = V e + V°, we end up with the expression for the interaction vertex: 

14^(k,k') = y {(iv v ) a / 3 (icr v )+ s + (i(k ■ cr)^)^^ • (t)ct v )+ s - (i[k x cr]a v ) a p{i^! x <r]cr y )+ 6 | . (S14) 

It has the form of Eq. (IS 6 I) with coefficients ao = a\ = 1, 02 = —1, 03 = 0. These are exactly the coefficients in the 
first row of Table EH 

We see that we have superconducting instabilities in s-wave and gyrotropic p-wave channels. As we mentioned 
above, they have identical transition temperatures, T c ~ woexp(—1 /zx|Vo|), where loq is some low-energy cutoff, and 
v is the density of states at Fermi level. The detailed derivation of the transition temperatures for the case of 
unconventional superconductors can be found, for example, in IS4l. lS5|. 

All the cases of other possible interactions in three and two dimensions can be analyzed absolutely analogously. 
The results of the decomposition into different superconducting channels are gathered in Table ED 


B. Multipolar interaction with reduced rotational symmetry 

Here we demonstrate that multipolar interaction in crystals with reduced rotational symmetry leads to the attraction 
in the correspondent particle-particle channel in both two and three dimensions. Moreover, in two dimensions, the 
critical temperature for the multipolar channels is the same as for s-wave channel. 

We consider the 3D case first. We perform the same analysis as in the previous section, but now we consider 
diagonal and off-diagonal parts of the multipolar interaction separately. Specifically, we consider Hamiltonians 


tfeff = 4a4'- q/3 «k' 7 a k _ q5 £(lf (2k - q)) Q ,(T«(2k' - q)W (S15) 

k,k',q i 

H eff = y E a ka a k'- q /3 ak '7 a k - q <5 - q))erf(lf (2k' - q))^, (S16) 

k,k', q ijtj 

where T^ are defined in Table El The whole analysis is similar to what we did before, with the only difference that 
now we need to consider diagonal and off-diagonal pairing channels separately: 

V afhS { k,k') = ^ (a o (ia*) a 0(*<7% + a 1 («Fi(k) ( 7») a/) (iF 1 (k')(7»)+ + a 2 (iF^k)a y ) afi (iF^k r )a v )+ 5 + 

"b ^3 ( k ) < ^)cd3(^-f 1 3^(k / )^"^) 7 (5 J ■ (S17) 

* i¥=j ) 

The results of the decomposition into different channels are presented in Table ism 

We see that, if considered separately, diagonal and off-diagonal parts of the multipolar interaction also lead to the 
attraction in the correspondent pairing channel. However, the interaction strength is smaller (though comparable) 
than that for the s-wave channel. As a result, the transition temperature for the s-wave channel is much higher. 

As we already mentioned, T^ is a traceless symmetric tensor, so it has only 5 independent components. That is 
why the diagonal part of T^ forms duplet (T|, T§) with 


r 


i 

3 


pxx x'W 
1 3 1 3 > 


r§ = 2 rf 


pii rvv 
1 3 1 3 > 


(S18) 


9 


type of interaction 

a 0 

Ol 

02 

at 

«3 

(Qa’WQs’t-q)) 11 

8/3 

0 

-2/3 

1 

-1 


4 

0 

-1 

-3/2 

1/2 


type of interaction 

a 0 

Ol 

02 

a 3 

a% 

a% 

0r(q)0r(-q) 

1 

0 

0 

-1 

1 

-1 

Qf(q)Qr(-q) 

1 

0 

0 

-1 

-1 

1 


TABLE SII: Decomposition of the multipolar interaction with reduced rotational symmetry into pairing channels. Diagonal 
and off-diagonal parts of the interaction in 3D are defined as {Q^ q)) d = q), (Q 3 7 (q)Q 3 *(^q))° = 

.. Q 3 7 '(q)Q 3 I ( — q) correspondingly, see Eqn. (IS15I) . AS16I) . 


rather than triplet (r 3 x , T| y , T 3 Z ). Hopefully, the diagonal part of the interaction can easily be rewritten in terms of 
(T 3 , r§) in both particle-hole and particle-particle channels: 

(r^(k)) a 4 (rr(k'))fl 7 + (rf (k)U(rf (k '))^ 7 + (rf(k)) a5 (rf (k'))^ = 

= 2(rl(k)) a5 (rl(k')) M + |(rl(k)) a5 (r2(k')W (si9) 


(*rr(k)^) a/ 3 (*rr(k >% + (*rf (k)^) a/3 (*rf (k')a^)+ + (ir 3 -(k)^)^(*rf (k')^)+ = 

= 2 (iri(kKM*ri(k>% + |(ir§(kK) a 0 (»r§(kV)+. (S 20 ) 

Finally, we consider the multipolar interaction with broken rotational symmetry in 2D. Unlike 3D case, T 4 splits 
now into two one-dimensional representations with form factor Tf^k) = k x a x — k y a v and T^k) = k x a y + k y a x . 
Consequently, diagonal and off-diagonal parts of interaction are given now by 

tfeff = a kc a k'- q ^k' 7 a k _ q5 (fr ( 2 k - q)W(fr( 2 k' - q))/ 3 7 > (S21) 

k,k',q 

H °S = T E «- q ^k' 7 a k - q5 (ff (2k - q)) a4 (ff (2k' - q))/ 3 7 - (S22) 

k , k ', q 

We emphasize that there is no any summation over different components of T 4 here. Indeed, due to the equalities 
ff^k) = — f^k) and f^k) = f| x (k), this summation would only lead to the overall factor 2 in front of the 
interaction Hamiltonian. Analogously, we do not double count different components when decompose the effective 
interaction into different pairing channels. 

The results of the decomposition are again gathered in Table ISIII In contrast to the 3D case, s-wave channel 
and corresponding multipolar channels now have the same interaction strength, and, as a result, the same critical 
temperature T c ~ u >0 exp(—l/V|Vo|). 

In order to obtain the last row of Table m one need to sum up diagonal and off-diagonal contributions from Table 
eh and to multiply it by 2. This factor 2 comes exactly from the double counting of the components of 1?4 that we 
mentioned above. 
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